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We propose a general procedure for reducing the three-dimensional Schrodinger equation for atoms 
moving along a strongly confining atomic waveguide to an effective one-dimensional equation. This 
procedure is applied to the case of a rotating closed-loop waveguide. The possibility of including 
mean-field atomic interactions is presented. Application of the general theory to characterize a new 
concept of atomic waveguide based on optical tweezers is finally discussed. 
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I. INTRODUCTION 

A very promising and challenging experiment to be 
performed in the near future using coherent matter waves 
is the observation of a rotation-induced supercurrent 
around a closed loop. This not only will attract a broad 
interest from the point of view of fundamental physics as 
a direct manifestation of superfluidity 0, but may also 
open the way to new kinds of high-precision rotation sen- 
sors based on matter- wave 2J , rather than light- wave 3j 
interferometry. The irrotationality constraint of super- 
fluids is in fact softened in the multiply connected geom- 
etry of a closed loop waveguide, which then appears as 
an ideal playground for the study of quantized vorticity 
and the related quantum interference phenomena. 

As it has been originally pointed out by Bloch , and 
then developed in more detail by Ueda and co-workers 
0, an appropriately tailored rotating potential can be 
used to transfer vorticity into a Bose-Einstein conden- 
sate trapped in a closed loop geometry. The 2D numeri- 
cal simulations reported in l5l| , which fully include the ef- 
fect of interactions, have demonstrated a close qualitative 
analogy between the physics of rotating condensates and 
the physics of condensates in ID optical lattices, an anal- 
ogy that can be made quantitative by correctly taking 
into account the modifications to the trapping potential 
due to rotation. Metastability appears for a rotating flux 
as a consequence of flux quantization: points of vanishing 
density have in fact to be introduced in the wavefunction 
if quantized vortices are to penetrate or exit the cloud 
0. When the velocity of the superfluid flow exceeds the 
velocity of sound, the system becomes energetically un- 
stable according to the Landau criterion of superfluidity. 
In this regime, the interaction of the flowing fluid with a 
stationary defect is able to create excitations in the fluid, 
which can eventually lead to dissipation of the supercur- 
rent by means of phase slippage processes 0, Q . 

In order to simplify the theoretical study of the gen- 



eration and dynamics of supercurrents in realistic loop 
configurations (e.g. a toroidal trap), it is very conve- 
nient to be able to isolate the longitudinal dynamics along 
the, possibly rotating, waveguide so to reduce the full 3D 
problem to ID. This is one of the central points of the 
present paper. 

Once the 3D problem is reduced to a ID one, the for- 
mal connection between annular rotating Bose-Einstein 
condensates and Bose-Einstein condensates in optical lat- 
tices becomes apparent, and one can start taking fully ad- 
vantage of the large amount of literature appeared on the 
latter subject. A remarkable example of this connection 
are the so-called swallow tails (instable parts of the band 
structure at the edge and in the centre of the Brillouin 
zone) shown by the band structure of a BEG flowing in 
an optical lattice 0, ■ ^ it has been pointed out in 
Ref. they play an important role in the nucleation 
of vortices and supercurrents in rotating BECs. Further 
very interesting phenomena that have been studied in the 
context of BECs in optical lattices are dynamical insta- 
bilities 10, 0, 0, and gap solitons 01 j which are 
expected to have interesting counterparts in the physics 
of rotating BECs. 

First observations ^3 of Bose-Einstein condensates 
trapped in toroidal traps have been recently reported in 
Refs. These experimental setups were based on 

magnetic potential, which introduces some limitations on 
the geometries that can be obtained and, more specifi- 
cally, on the flexibility of the setup with respect to time 
and space modulations of the confining potential. 

The other central point of our paper is in fact to pro- 
pose a novel realization of toroidal atomic waveguide, 
that makes use of the attractive optical potential of a 
red-detuned laser beam as an optical tweezer The 
main advantage of our proposal with respect to recent 
related ones 21] is the flexibility of the setup. A toroidal 
trapping potential with a strong transverse confinement 
is obtained by rapidly moving the focus point of the laser. 
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As the trajectory of the focus point can be arbitrarily 
chosen, as well as its speed of motion, any shape can 
in principle be obtained for the ID waveguide, and any 
longitudinal potential can be applied onto the atoms in 
addition to the transverse confinement. Furthermore, the 
shape of the curve can be changed in time, so to obtain, 
e.g., rotating waveguide traps. 

The paper is organized as follows. In Sec. ^ we put 
the problem on a precise mathematical basis and then 
in Sec, mil we introduce our decoupling scheme to reduce 
the three-dimensional Schrodingcr equation to an effec- 
tive one-dimensional one. Sec. lIII Al and Scc. lIII Bl discuss 
the effect of a longitudinal dependence of the transverse 
trapping frequency, and the effect of the interatomic in- 
teractions. The main theoretical results of this paper are 
given in Sec. II VI where the decoupling scheme is gener- 
alized to the case of a rotating waveguide. In Sec. 
the theoretical approach is applied to our proposal of ID 
waveguide based on a rapidly moving optical tweezer. 
Conclusions arc finally drawn in Sec. IVII 



II. 



QUANTITATIVE DESCRIPTION OF THE 
ATOMIC WAVEGUIDE 



Consider an atomic waveguide whose axis follows a reg- 
ular curve C parametrically defined by the vector rc(s), 
s being the arclength coordinate along C. At each point 
of the curve, we can define the Frenet frame (t, n, b) as 
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where n and r are respectively known as the curvature 
and the torsion of C In the vicinity of C, a local 

system {s,u,v) of coordinates can be introduced, such 
that 



r(s, u, v) = rc(s) -I- uN(s) -I- uB(s). 



(5) 



The transverse frame (t,N, B) is related to the Frenet 
frame (t, n, b) by a simple rotation around t 



cos ti 
— sin( 
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With this choice of coordinates, the gradient has the fol- 
lowing simple form 



V = th-^ds + 'Ndu + Bdy 
where the scale factor 



h{s , u, v) = I — K (u cos ( 



V sm f 



(7) 



(8) 



depends on the torsion r only via the angle 9. 

The transverse confinement in the (N, B) plane orthog- 
onal to the waveguide axis is assumed to be harmonic and 
of the form 



V±{U,V) = y (ujIu^+UjIv^) , 



(9) 



where M is the atomic mass and the transverse 

trapping frequencies in respectively the N and B direc- 
tions (which can depend on the longitudinal coordinate 
s). For the sake of simplicity, the discussion that fol- 
lows will be restricted to two mostly significant cases. In 
Sec. mil the curve is allowed to have a non- vanishing tor- 
sion T, but the transverse trapping is taken as isotropic 
This condition is enough to rule out ef- 
fects due to the finite torsion 0. In the Sees II VI andlVia 
different situation is considered, where the curve is taken 
to belong to the plane orthogonal to the rotation axis, 
while the trapping potential can have different frequen- 
cies 7^ 'jJv in the two orthogonal directions respectively 
in and perpendicular to the plane. 

The key assumption of our treatment is the strong con- 
finement hypothesis, where the extension 



Muji_ 



(10) 



of the transverse ground state is much smaller than all 
typical length scales of the curve C, namely 



Kf7 < 1 



\n \o- < k 



k'V« 



(11) 



Here, and throughout the whole paper, primed quantities 
denote derivation with respect to the longitudinal coor- 
dinate s. The reason why conditions Hll|l involve up to 
the second derivatives of the curvature will become clear 
in the following. These conditions also guarantee that 
the coordinate system (O can be safely used to describe 
the transverse extension of the wavefunction. 

If the waveguide is at rest in a reference frame rotating 
at an angular speed Jl, the stationary Schrodinge r equa- 
tion in the rotating waveguide has then the form 



of an angle 9 such that 



A9 



(6) 



- - f2 • (r X p) 



/i^(s, u, v) = 

+ V1(m,w) +yoxt(s) 



^'(s, u, v) 



(12) 
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where the momentum operator is defined as usual as p = 
~ih\J . In the (s,m, w) coordinates, the gradient has the 
form 0. The term T4xt(s) describes any weak potential 
acting on the atoms along the direction of the waveguide 
in addition to the waveguide trapping. 

Note that the normalization of the wavefunction 'J has 
to take into account the new metric associated to the 
{s,u,v) coordinates 

ds du Av /i(s,u,i;) |*(s,it,w)p = 1 . (13) 

As h{s, u, v), defined in is not factorizable as a prod- 
uct of functions of respectively the longitudinal s and 
transverse (u, v) coordinates, this condition is not conve- 
nient for decoupling the transverse and the longitudinal 
dynamics. It is then useful to introduce a rescaled wave- 
function $ = whose normalization is the usual one 

dsdudu |<i>(s, u, w)p = 1 . (14) 

The Hamiltonian operator acting on the wavefunction 
$(s, u, v) then has the form 



freedom are eliminated and summarized as an effective 
potential acting on the nuclei. 

For any longitudinal wavefunction g{s), we define a 
transverse Hamiltonian h at the position s such that 

hfs = Hifsg). (17) 

The key assumption of our approach is to assume the 
transverse motion to be frozen in the ground state of 
h. Our aim is to reduce the three-dimensional problem to 
a one-dimensional one, by integrating over the transverse 
degrees of freedom 



W(s) = jj dudv,C-{u,v)Hf^iu,v)g{s) = Hgg{s), 

(18) 

hence eliminating adiabatically the transverse motion. 
An explicit form of Hg can be obtained by means of 
a perturbative expansion by separating in h the differ- 
ent contributions due to the transverse and longitudinal 
degrees of freedom, namely 





H = -- 



2M 



ds{h-'is,u,v)ds) + dl^ + dl 

5[/i'(s,u,t>)]2 h"{s,u,v) 



4/i^(s,u, t;) 4:h'^{s,u,v) 2h^(s,u,v) 



(15) 



— \/h{s,u,v) fl • (r X p) 
+ Va_{u,v) + Vc^t{s) . 



In the next sections, we shall proceed with the decoupling 
of the transverse and the longitudinal dynamics under the 
strong confinement hypothesis. 



III. DECOUPLING PROCEDURE IN THE 
NON-ROTATING CASE 

In this section, we shall start by considering the sim- 
plest case of a non-rotating waveguide {fl = 0) with 
a spatially constant trapping frequency {uj'j_ =0). It 
is useful to rewrite the wavefunction $ as the product 
^{s,u^v) = fs{u^v) g{s) of a longitudinal wavefunction 
g{s) and a transverse wavefunction fs{u,v) in general 
dependent on the longitudinal coordinate s. The nor- 
malization conditions can then be written as 



h = HQ + W 



with 



Ho = g{s) 



2M 



2 

uu 



-dt,)+V^{u,v) 



and 



W = 



'2M 



h-\s,u,v)[g{s)dl + 2g'{s)ds 



g{sW{s) ^ 5g{s)[h'{s,u,v)r 



4/i^(s, u, v) 
g{s)h"{s,u,v) 



2h^{s, u, v) 



f g(s)Kxt(s) 



(19) 



(20) 



(21) 



We treat perturbatively the Hamiltonian W with respect 
to Hq, which a part from the multiplying factor g, cor- 
responds to an harmonic oscillator with energy /lo^x. All 
terms of W are much smaller than |(7|fitij^ provided 



dudv\Mu,v)\ =1, \g"\«\g\/<y' and \g'\«\g\/a, (22) 

r (16) 

/ ds|(7(s)| — 1. conditions which will be self-consistently verified at the 

end of the procedure, thanks to conditions (|ll|l . 
We now proceed in the spirit of the Born-Oppenheimer Perturbation theory at first-order (in the small pa- 
approximation 23] where the fast electronic degrees of rameter k^ct^) allows to replace /° in Eq. H18|) with the 
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ground state wavefunction /o of the harmonic oscillator 
of frequency wj^, given by 



20 



VTTcr 



(23) 



leading to the following effective one-dimensional 
Schrodinger equation for the longitudinal wavefunction 
9{s) 



'2M ds2 



8M 



9- (24) 



We easily recognize the usual kinetic energy term, the 
zero-point energy of the two-dimensional transverse trap- 
ping potential, and the weak external potential Vl^xt- The 
last term of equation 1)24(1 gives an effective potential 
proportional to the square of the curvature as discussed 
in 13 ■ The smoothness of the waveguide, as quantified 
by conditions lfTT|l . guarantees that conditions are 
satisfied and hence the self-consistency of the approach. 

This result can be illustrated on the specific example of 
an elliptical waveguide, whose axis are respectively equal 
to i? cosh 770 and i?sinh77o. The parameter 770 character- 
izes its eccentricity (the larger 770, the closer to a circle) 
and R gives the overall scale. The strong confinement 
hypothesis is then a <ti -Rsinh?7o. The one-dimension 
equation is 



W(s) 



- fliO 



2Mds2 ' " ^ 8M 
The curvature k has the simple expression 



g{s) . (25) 
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Figure 1: Result of the numerical integration of the 2D 
Schrodinger equation in the case of an elliptical waveguide 
with strong transverse confinement. Parameters: semiaxis 
b/a ~ 34.5; a/a ~ 15.6. 




^ ^ coshr/psinhr/o ^ ^^^^ 
i?(sinh^ rjQ + sin^ w)^/^ 

where w is the parametric angle along the ellipse, related 
to the arclength coordinate s by the differential relation 



ds = i?Y sinh^ rjQ + sin^ w dw . (27) 

The curvature k has its maxima on the great axis of 
the ellipse. The curvature-induced effective potential is 
thus minimum at these points. Its effect is illustrated 
in Fig^ where the results of a numerical integration of 
the 2D Schrodinger equation in cartesian coordinates are 
shown (the third dimension was neglected for the sake 
of simplicity). Due to the curvature- induced potential, 
the atomic density is maximum on the great axis of the 
ellipse. 

A more quantitative comparison between the full 2D 
Schrodinger equation in cartesian coordinates and the 
effective ID Schrodinger equation H24|l is shown in Fig[21 
for the case of an elliptical waveguide with strong con- 
finement. The agreement is remarkable. 



Figure 2: Result of the numerical integration of the 
Schrodinger equation in the case of an elliptical waveguide 
with strong transverse confinement. The graph shows the ra- 
tio between the density at the extrema of the two orthogonal 
axes of the ellipse, as a function of the ratio between the length 
of theses axes b/a. The points come from the simulation of 
the 2D Schrodinger equation in cartesian coordinates, while 
the full line comes from the ID simulation of equation 12411 . 



A. Effect of a longitudinal variation of uj± 

The case of a transverse trapping frequency lo±{s) with 
a non-trivial dependence on the longitudinal coordinate 
s is addressed in the present section. This induces a 
longitudinal variation of the transverse wavefunction /os 
and introduce new terms in Hg due to the non- vanishing 
longitudinal derivatives of fos ■ Applying the same proce- 
dure as in the previous section, and limiting ourselves to 
the first order in k^ct^, we get the following ID effective 
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Figure 3: (a) Result of the numerical integration of the 2D 
Schrodinger equation in the case of an straight waveguide with 
a strong confinement of spatially- dependent transverse con- 
finement frequency u>±{x) — Lu±y^l + 4:{Suj±/ll!±) cos{ttx/\). 
The colour code is the same as in FigQ Parameters: 
Suj±/ijj± = 0.1, X/a ~ 22.4. An external potential such that 
Vcxt + Vl — MLj±{x)'^y^ /2 — huj±{x)/2 has been introduced in 
order to compensate the spatial modulation of the zero-point 
trapping energy, (b) comparison between the 2D density in- 
tegrated along direction y (red dots) and the solution if the 
effective ID Schrodinger equation (full line). 



equation for g(s) 



= 



2M ds2 



8M IGMojj 



(28) 

The dependence of uj± on s not only gives an s-dependent 
potential energy equal to zero-point energy hu!±{s), but 
also adds a further contribution proportional to w^/w^. 
Consistency with our decoupling procedure requires that 
^ Kuj± and jw" | ^ k^uj^. This implies in particular 
that the term proportional to o^^/lj^ is much smaller 
than the zero-point energy term. 

In order to check the effect of the new potential term 
proportional to w^/w^, we have numerically solved the 
2D Schrodinger equation in cartesian coordinates for the 
case of a straight waveguide with a modulated uj±{s). In 
order to put the effect of the uj'^/ui'j^ in better evidence, 
the external potential V^xt has been chosen in such a way 
to compensate the modulation of the zero-point energy of 

the transverse harmonic trapping /ia;_L(s)/2-|-14xt(s) — 0. 
As shown in FiglSJ the density modulation has indeed 
the same periodicity of uj'^ / and is in quantitatively 
good agreement with the numerical solution of the one- 
dimensional equation 



B. Effect of interatomic interactions 

In a 3D geometry, mean-field interactions are included 
in the the Gross-Pitaevskii equation Bl as nonlinear 
terms of the form: 



Vint = gaoNa 



fs{u,v)g{s) 



(29) 



being g^u = Anfi^asu /M, a^D the three-dimensional 
scattering length and A^o is the total number of atoms 01 ■ 
The Gross-Pitaevskii equation in presence of a periodic 
potential has been extensively studied in the context of 
Bose-Einstein condensates in optical lattices. The issue 
of the factorization of the wavefunction in its transverse 
and longitudinal part in the presence of interactions is in 
general a non trivial one. 

In the language of the present paper, the presence 
of interactions requires that the transverse wavefunc- 
tion /o(m, w) is now solution of the ground state Gross- 
Pitaevskii equation 



2m 



{dL+dl) + Nog3D\g{s)\'\fos{u,v)\^+ 

+ V±(u,v) /os (u, w) = M_L fos {u,v) (30) 



which results from the inclusion of the mean-field inter- 
action term into the transverse Hamiltonian (|2()|l . Here 
fi± has the meaning of a transverse chemical poten- 
tial. In the general case, Eq. H30|l has to be solved self- 
consistently with the equation determining the longitu- 
dinal wavefunction g{s), which can be a computationally 
intensive task unless clever schemes are adopted, as e.g. 
in 24]. 

A very simple formulation can be obtained in the lim- 
iting case of very strong radial confinement, when in- 
teractions have a negligible effect on the shape of the 
transverse ground state wavefunction /o(u, u). In this 
limit, interactions only provide a mean- field energy term 
proportional to the local density to be included in the 
longitudinal equation. This has then the usual form of 
a ID Gross-Pitaevskii equation including the curvature 
term discussed in the previous section 



2M ds2 



+ K.xt(s)(7(s) 



8M 



g{s) + giDNo\g{s)\^gis), (31) 



where the effective ID coupling constant is defined as: 



Ma2 



IV. CASE OF A ROTATING WAVEGUIDE 

We now turn to the more general case of a rotating 
waveguide. For simplicity, we shall from now on assume 
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that the curve C is included in the plane {x,y), that the 
rotation vector fl is orthogonal to this plane, and that 
the origin of rc coincides with the centre of rotation. We 
shall furthermore consider in what follows the case of a 
constant uj±_ independent of the position on the curve. 

Repeating the same steps as in Sec. lIIII the transverse 
Hamiltonian can be decomposed as 



h = Ho + W 



(33) 



with 



2M ' 



- TTTJ^vv + gVi_ , (34) 



and 



W = W + iMl 
_ grnh' 



2h 



gu 



guh! 



rng 



--MQ^^g. (35) 



In the last expression, the following shorthand notations 
have been used: rt — rc ■ t and r„ = rc • n. The ground 
state of Hq is now given by 



where we have used the identity dr„/ds = Krt- As com- 
pared to the non- rotating case of H24|l . additional terms 
appear here due to the rotation. The first one is the 
one-dimensional form of the gauge term appearing in the 
kinetic energy term in the rotating frame and mostly af- 
fects the phase of the wavefunction. The second rotation- 
induced term is the classical centrifugal energy term. 
This term, along with the curvature-induced term, the 
external potential term, and the zero-point energy term, 
can be used to transfer angular momentum from the trap 
to the atoms and then establish a finite phase circulation 
in the condensate. 



A. Analogy with optical lattices 

Consider for simplicity a circular waveguide of radius R 
in the presence of a rotating periodic potential of period 
2tt/£ {i is integer) in the angular coordinate 6 



v{e,t) = Vocos[e{e -nt)]. 



(40) 



Once mean-field interactions are included in the same 
way as done in 1)31(1 . equation 1)39(1 can be casted in the 
simple form 



(36) 



Assuming a moderate rotation speed <^ lu±_ and 
MfJ^r^ <C fiuj±, an explicit calculation of the longitudi- 
nal derivatives of fos shows that the following inequalities 
are satisfied 



5/, 



Os 



ds 



fos 

a 



and 



d'hs 



fos 



(37) 



which guarantee that W can be treated as a small per- 
turbation with respect to Hq. To the first order in K^a'^, 
one then has 



dudvf;^iHo + W)fos 



2M- 



8M 



ft 

2 * 
9 + rng' 



(38) 



and the ID equation for the rotating waveguide can fi- 
nally be written as 



1 

2M 



ih h Af t • (rc X n) 

ds 



2^2 



hiui 



8M 



1 



2„2 



5, (39) 



80 



Vocos{i9)+g\ip\' 



^{6) = ~M0) 



(41) 

which shows a strong formal analogy with the problem 
of a Bose-Einstein condensate in a ID optical lattice 
with periodic boundary conditions |^. Here, ip{0) is 
the longitudinal wavefunction, g is the coupling constant 
due to interatomic interactions, fl is the chemical po- 
tential shifted by constant potential terms and all ener- 
gies have been expressed in units of fi^ /{2MR^). The 
quantity fi/2 plays the role of the quasi-momentum for 
a Bloch wave in the periodic potential of the lattice, 
with a subtle difference arising from the different peri- 
odic boundary conditions: in the periodic potential of 
the lattice, the allowed values of quasi-momentum are 
integer multiples of h2-K/Nd, d being the lattice spac- 
ing and N being the number of lattice wells present in 
the system. On the other hand, all values of are al- 
lowed in the present case of a rotating waveguide, the 
single- valuedness of the wavefunction ip[9) being ensured 
by the Bloch theorem on the whole length of the ring, 
i.e. Lp{9 + 2-^/1) = eyi^\i2mT / P\tf{9) with n integer. The 
complete band structure is generated by the set of values 
n ~ 0,1, 2,...^ — 1, giving rise to (, independent sub- 
bands periodic in with period 21. 

The circulation of the different states is a function of n 
and VL (for instance the lowest energy states at 2ri — 1 < 
Q. <2n-\-l have circulation n). The ^-fold modulation of 
the potential along the waveguide allows only mixing of 
states at circulation n with states at circulation n ± £, in 
correspondence of the sub-band gaps. Hence, following 
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adiabatically a given sub-band by very slowly increasing 
(or decreasing) the rotation frequency, one can transform 
a state at a given circulation n into a state at circulation 
n±(.. 



V. EXPERIMENTAL ISSUES 

A possible way of achieving experimentally an annu- 
lar condensate with strong transverse confinement is to 
use a magnetic toroidal trap, as reported in 18, 19]. In 
this section, a completely different experimental configu- 
ration is proposed, and its advantages pointed out. For 
the sake of simplicity, we shall concentrate our attention 
to the most relevant case of a planar waveguide whose 
axis lays on the horizontal xy plane. A strong confine- 
ment in the vertical z direction can be obtained using a 
horizontal light sheet which provides a potential Vz{z), 
while the curvilinear waveguide profile can be created 
by rapidly moving the focus point of a red-detuned laser 
beam (acting as an optical tweezer 0|) along the curve 
C, at a, possibly time-dependent speed v{t). This can be 
obtained e.g. by reflecting the laser light onto vibrating 
mirrors or using acousto-optic modulators. If the move- 
ment of the focus point is much faster than the transverse 
trapping frequency, then the atoms will see the following 
averaged potential 



f dt 

V{x,y,z) = V,{z) + j;VtUx-xc{t),y~ycit)), (42) 

where Vz{z) is the trapping potential due to the light 
sheet and Vtw{x, y) the trapping potential due to the op- 
tical tweezer. The pair {xc{t),yc{t)) defines the position 
of the laser focus on the xy plane at time t, which spans 
the curve C in the period T {x and y are here cartesian 
coordinates). The tweezer potential is attractive, and 
can be taken as having a Gaussian shape of waist w and 
peak value < 



Ftw(x,y) = C exp[-{x'+y^)/w% 



(43) 



It is then easy to obtain expressions for the potential 
terms V±_ and t4xt appearing in the Hamiltonian (|12|l 



Kxt(s) 



UJ. 



Tv{s) 
Tv{s) 



with wH, 



,(44) 
(45) 



A remarkable fact has to be noted: both w^(s) and 
Vcxt (s) are inversely proportional to the speed of the laser 
focus when passing in the neighbourhood of the point s 
under examination, and to its overall orbital period. A 
possible way of adding a spatially modulated external po- 
tential along the waveguide is therefore to simply modu- 
late the speed of motion of the focus along the curve. 



Plugging in (|44|l typical values of intensity / = 100 W 
and detuning S = 300 nm for optical tweezers with 
w = 30 fim, one can see that transverse extensions as 
low as 1.16 /im could be achieved on a 150 /im radius 
torus with such a device, ensuring the validity of the 
strong confinement hypothesis. The main advantage of 
this method over the magnetic traps studied e.g. in T^l 
is that any shape of the curve C can be obtained without 
any additional difficulty, and this can also be modified 
in the course of the experiment in order to obtain e.g. a 
rotating waveguide. 

Another important requirement for the study of super- 
current to be possible is that the longitudinal potentials 
are not strong enough to fragment the condensate. As 
an example, it is interesting to estimate the effect of a 
tilting of the light sheet that is used to vertically confine 
the atoms. For a tilting angle a from the horizontal, the 
potential energy difference at the extrema of the waveg- 
uide due to gravity is equal to 2MgQRsma, where go is 
the gravity field acceleration and 2R is the horizontal dis- 
tance between opposite points of the waveguide. Within 
the Thomas-Fermi approximation, fragmentation occurs 
if the gravitational energy difference is larger than the 
mean-field interaction energy, that is if 



2MgGRsma > 



2ttR ■ 



(46) 



For a system of Nq = 10^ atoms of ^^Rb with scattering 
length a — 5.77 nm in a waveguide of radius R = 150 fiia 
and a = 1.16 /im, the system remain connected provided 
a < 0.15°, which is a rather stringent but not unreach- 
able experimental requirement. Note that the gravita- 
tional potential does not rotate with the waveguide when 
this is set into rotation, so that it might possibly act as 
a defect moving through the one-dimcnsional fluid 0, 0] . 

In the case of very tight confinement and frozen trans- 
verse dynamics the issue of phase fluctuations in the ID 
condensate comes into play. We believe that the rota- 
tional properties of the condensate should not however be 
disturbed at least in the meanfield regime at low enough 
temperature. 

Another issue of great importance from the experi- 
mental point of view is the possibility of measuring the 
number of vorticity quanta present around the annular 
condensate. A measurement after expansion has been re- 
cently predicted to be able to provide a clear answer [23| , 
but a non-destructive measurement would be preferable 
in view of applications as a rotation sensor. In the case of 
a rotating non-circular (e.g., elliptical) waveguide, the de- 
formation of the periodic density modulation due to the 
centrifugal potential should in principle provide a way of 
measuring the supercurrent. If this signal is too weak 
to be detected, one could still resort to other techniques, 
e.g. the analysis of collective modes or the measure- 
ment of the momentum distribution by means of Bragg 
spectroscopy j29] or slow-light imaging [30l |. 
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VI. CONCLUSION 

In this paper, we have developed a formahsm which 
is able to reduce the three-dimensional problem of 
atomic propagation along an atomic waveguide to one- 
dimensional equations. Under a strong transverse con- 
finement hypothesis, the transverse extension of the 
wavefunction is much smaller than the curvature radius 
of the waveguide and the wavefunction can be factorized 
into the product of its longitudinal and transverse parts. 
Our formalism is then applied to a novel concept of opti- 
cal waveguide which combines the possibility of having a 
strong confinement with a great flexibility in the design 
of the, possibly rotating, waveguide shape. 

Such a description provides a simple yet accurate start- 
ing point for analytical studies and numerical simula- 
tions, as well as for the design of experimental setups. 
Our framework is in fact able to considerably simplify 



the theoretical analysis while still keeping track of the 
relevant degrees of freedom in a quantitative way. Possi- 
ble applications range from the determination of the best 
experimental sequence to nucleate a supercurrent along 
a ring-shaped Bose-Einstein condensate, to the study of 
the response to a global rotation in a sort of matter-wave 
gyroscope. 
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